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Abstract 

We calculate the spectrum of the relic gravitational wave due to the trans-Planckian effect in 
which the standard linear dispersion relations may be modified. Of the modified dispersion relations 
suggested in literatures which have investigated the trans-Planckian effect, we especially use the 
Corley-Jacobson dispersion relations. The Corley-Jacobson type modified dispersion relations can 
be obtained from Hofava-Lifshitz gravity which is non-relativistic and UV complete. Although it is 
not clear how the transitions from Hoi^ava-Lifshitz gravity in the UV regime to Einstein gravity in 
the IR limit occur, we assume Hofava-Lifshitz gravity regime is followed by the inflationary phase 
in Einstein gravity. 
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I. INTRODUCTION 



In spite of the success of the inflationary scenario, inflation is still encountered several 
unsolved problems. Among them, the trans-Planckian problem [l| is related to the funda- 
mental problem of the high energy physics i.e quantum theory of gravity. Unfortunately 
there is no successful quantum gravity theory for handling the spacetime around the Planck 
scale, so one way out of such difficulties is to use the effective modified dispersion relations to 
mimic the effect of the trans-Planckian physics. Of the several modified dispersion relations 
suggested in many literatures, we will use the Corley-Jacobson dispersion relation m 
this work. 

The Corley-Jacobson type dispersion relations could be derived naturally from the re- 
cently proposed quantum gravity model called Hofava-Lifshitz(HL) gravity sj 4]. HL gravity 
is a non-relativistic and renormalizable theory in the UV region. In addition, because of the 
different scaling behaviors of the time and space (t I'^t, x /x), while the kinetic term is 
still quadratic order in time derivatives of the metric in the action, the potential terms may 



have up to 6th order spatial derivatives [5|_[6 
the large scale structure without inflation 



HL gravity may also provide the seeds for 
7|. On the other hand, it is still controversial 
if HL gravity is a complete theory[8|] which would be dependent upon the projectability and 
the detail balance condition sl. No explicit inconsistencies has been pointed out against the 
version without the detailed balance condition with the projectability condition |9|. And it 
is also not clear how the transition from HL gravity in the UV regime to Einstein gravity in 
the IR limit takes place. 

In this paper, we calculate the spectrum of the relic gravitational wave in order to in- 
vestigate the effect of the trans-Planckian physics. Especially, we use the Corley-Jacobson 
dispersion relations taking HL gravity into account. For this purpose, we assume the HL 
gravity phase is followed by the inflationary phase in Einstein gravity. Though the transition 
time from the trans-Planckian regime to the usual standard inflationary period is usually de- 
termined by the cut-off momentum, we will leave it as a free parameter assuming it keep the 
information about the transition mechanism. The resulting spectrum could be constrained 
from the experimental bounds such as the ground based gravitational wave detector (LIGO, 
VIRGO) which covers the frequency range lOHz < / < lO^Hz and future planned space 
based one (LISA) which is expected to cover 10~^Hz < / < IHz or the cosmological bounds 



(nucleosynthesis, CMB) in the low frequency range 10|. 

The paper is organized as follows. In Sect. [TTl we calculate the power spectrum of the 
tensor perturbations with the Corley-Jacobson dispersion relations. In Sect. IIIIl we obtain 
the equation for the tensor perturbations in HL gravity and calculate the power spectrum. 
In Sect. llVt we investigate the effect of the trans-Planckian physics and HL gravity on the 
spectrum of the relic gravitational wave. And finally, we summarize in Sect. |Vl 

II. GRAVITATIONAL WAVE SPECTRUM IN THE TRANS-PLANCKIAN 
PHYSICS 

We briefly review the effect of the trans-Planckian physics on the power spectrum in this 
section. With the metric 

rfs^ = a\r])[-d7]^ + {S,j + h,j)dx'dx^], (1) 
where hij is a transverse, traceless perturbation which can be expanded as 

^''^""^ = 11 (^F^ S ^^Mv)e'''■'', (2) 

and eij is a polarization tensor, the tensor perturbations are governed by the equations 

ul+ -^^Uk = 0, (3) 

where a prime denotes a derivative with respect to the conformal time and Uk satisfy the 
normalization condition 

UkUk - u^u^ = I. (4) 

This equation is believed to hold only at the low energy scale, kp^ys -C Mpi, where kphys = 
k/a. As we go back in time around the Planck scale, quantum effects would be important and 
it seems to be necessary introducing the quantum gravity theory in order to describe correctly 
the Planck scale phenomena. Unfortunately, however, no successful quantum gravity theory 
exist at present. Though the string theory is believed to be a strong candidate of the 
quantum gravity theory, it is not complete yet. Thus, one way to overcome these difficulties 
and to mimic the quantum gravity effects around the Planck scale is to consider the modified 
dispersion relations. 



Taking into account of the modified dispersion relations, uj{k/a), the equation for Uk 
could have the following form 



a" 



ul+\^\kla)--juk = Q. (5) 

In this paper, we take the Corley- Jacobson dispersion relations 21] 11] as an example of the 
modified dispersion relations 

aM ^ 

where M is a cutoff scale and hm is an arbitrary coefficient. Here /c is a comoving wave 
number. Depending on the sign of hm which can be positive or negative we may reach 
to the different results. If < 0, uj"^ becomes a complex function for k > ]6m]~^/^'"aM and 
the WKB approximation may not be valid llj] in this region. For the purpose of this work. 



u\k/a) = k'(l + bm{47] )' (6) 



bearing in mind HL gravity, it is enoug. 
consider bm > 0. (For bm. < 0, see Ref. 
If we assume exponential inflation 



1 to limit to the positive bm- From now on we only 
ll)- 



the equation (jS]) yields 



where 



< + (e + <T^k)e'"*Y'" - ^ I ft = 0, (8) 



(rr \ 2m 
m) ■ 

For k > aM, since uj'^{k/a) ~ cr^/c^"*"'"^?]^'", the equation ([8]) becomes 

<+(aV"+V™-^)^fc = 0, (10) 
and its solutions are given by 



'/'Hl'\z) + c,J-—^Av\'^'Hi^\^), (11) 



'''^/ 4(1 + 771)'"^' ' ^^V4(l + m) 

where h!)^^ and hI^^ are the Hankel function of the first and second kind with the order z/, 
respectively, and 

-ikvr"", ^ = 7T77^- (12) 



1 + m' " 2(l + m) 



For bm < 0, the solutions are given by the modified Bessel functions. 

If we choose the positive frequency modes as initial conditions when k > aM, we can 
choose Ci = 0, C2 = 1 taking the normalization condition (jl]) into account. This choice of 
the initial vacuum is equivalent to the minimizing the energy density |l| up to the phase 
factor 

MVi) = u'kivi) = -^a/^' (13) 



if the adiabatic approximation {\uj'\ ^ u^) is hold. Then the mode solutions of ( ITOl) for 
k > aM can be written as 



\v\'/'Hi^\z). (14) 



4(1 + m) 

For k < aM, the dispersion relations recover the standard linear relations, a;^ ~ /c^, as in 
([3]) and then the mode solutions are given by 



Uk{r]) = ^ 1- 6-*'='' + -^ 1 + — e*'=^ (15) 



where the coefficients satisfy the normalization condition 

\Ak\'-\Bk\' = l. (16) 

The coefficients are determined through the matching conditions at ?7 = ?7i where t]i is 
the time when the modified dispersion relations take the standard linear form, 

, , , 1 ^ 

Through the matching conditions at ?7 = r/i, we obtain the coefficients 



where 



With these coefficients, we can compute the power spectrum of the gravitational wave at 
rj or krj <^ 1 

k^ M 

= (^^^ (^1 + 2\B,\' - Re{A,Bl)y (21) 

where we have used = at the second hne. 

For k\rii\ ^ 1, using f|T7j) and fl20l) . we finally obtain 

V,{k) ^ (^)(2-^|sin(2fc|,,|)), (22) 

where we have used the asymptotic form of the Hankel function 

/Ae-(-(-+|)|), if^»l, (23) 
V vrz 

and 

d V 

dz z 

The power spectrum ( l22l) is the well-known result of the trans-Planckian effect. The 
leading correction term due to the modified dispersion relations leads to ^. The detectability 
of this effect is an important concern at the future experiments. 



III. GRAVITATIONAL WAVE IN HORAVA-LIFSHITZ GRAVITY 



The quadratic action of Horava-Lifshitz gravity with the metric ([T]) is 

doSr = I dtd^xo? 



K 



(24) 



a-" " a" 

where Vi is a covariant derivative with respect to Sij and A = The coefficients in 

the action corresponding to the original theory ^ are 



ai 



"3 = TT^, "4 



8(1 -3A)^" 8 ' 2u;2' ' 2w^' ^^^^ 

where k, A, w and n are coupling constants with scaling dimensions [k\ = [w] = [A] = 0, [/i] = 
1 and A is the 3-dimensional cosmological constant with scaling dimension [A] = 2. Einstein 
gravity is recovered if A is approaching unity in the IR limit. 



The background Friedmann equation in HL gravity [3[ is given by 



where p is the background energy density of the matter and 



2, ,2 



cr 



3A^ (27) 



8(1 -3A) 

If p = 0, there is no cosmological solution In order to get de Sitter solutions, we 

need to introduce the matter where we assume the vacuum energy density of the matter is 
dominant. In what follows, we assume po = > and then the scale factor is given by 

^ (28) 



Hrj 



Note that by the analytic continuation of the parameters p — > ip and w"^ iw"^ 13| . 



12|. The sign of 



the de Sitter solutions can be obtained without introducing the matter 
the coefficients in the transformed action is also changed. In order to couple matter with 

n 

gravity consistently, the relativistic formulation of HL gravity are also investigated p^, in 
which, although being fully-relativistic, it results to be locally anisotropic in the time-like 
and space-like directions defined by a family of irrotational observers. 

In the action ( l24l) . because of term the tensor perturbations are not invariant under 
the parity transformation (x — * — x) [iSj. The chirality of the gravitational wave seems 
to give rise to the important effect for discriminating HL gravity from the other quantum 
gravity theories at the future experiment |l^[l6|. But since the term is subdominant 
compared to the term in the UV limit, we will drop such term in the present work. 

Varying the action fl2^ leads to the equation for the tensor perturbations in HL gravity 



u'l 



Uk = 0, (29) 



where we have used the mode expansion ([2]). In the UV limit, the dispersion relation of this 
equation is exactly equivalent to the Corley-Jacobson dispersion relation if m = 2. If only 
the dominant term is kept with putting 



= ^— , = 1, (30) 



then the positive mode solution is 



2 



«^ = \lY^\vrH^iM^)^ z = j {kr^r. (31) 



Since the Hankel function with the order z/ = 1/2 can be expressed as 




2 

TXZ 



= ^\li:e~'^ (32) 



the power spectrum of the tensor perturbations in HL gravity 

Vuik) <x k^lukl"^ <x (33) 

can be scale invariant in the UV region. If HL gravity can provide the seeds for the large 
scale structure without the inflationary phase as was suggested in , it would also be 
interesting to find some observational signature due to the tensor perturbations in HL gravity 
without introducing the inflationary phase But in this work, we assume HL gravity in 
the UV region is followed by the inflationary phase in Einstein gravity. 

We, then, move to the IR limit. Since it is expected Einstein gravity would be recovered 
in the IR limit with A = 1, the equation for the tensor perturbations takes the standard 
form as like in ([3]). The solutions, then, are 

u, = 4^(l-^) e~^'^ + + e^^\ (34) 

v^V kiqj V2A?V kiqj 



As in Sect. [TTl we use the matching conditions at 77 = r/i in order to determine the coefficients 
Ak and Bk- Although it is unclear how the transition from HL gravity in the UV region to 
Einstein gravity in the IR limit occurs, we assume that the transition occurs instantaneously 
at ?7 = 771. 

We, then, obtain the coefficients via the matching conditions 

A^;iw.,i)-(i.^-^)e-«<^^^^^ 

^1 + T^V""4j?!?i(--i). (35) 





ik\r]i\y/^\ k\r]i\J dz V 



e-^'^^'TT^I^l^i), (36) 



where 



Since, for kr] 1, 



2 

cos(2A;|?7i|) + — — - sin(2A;|?7i) 



/c^|?7iP k\r]i 
the power spectrum of the gravitational wave in the late time is, from (l2T!l . 



(39) 



where, similarly as in (fT7j) . we have used 

Although we use the relation (14T|) in order to get ( l40l) . since we do not have any informa- 
tion about the transition mechanism, we will, for the time being, leave rji as a free parameter 
which is assumed to keep the information about the transition. 

Then, for A;|?7i| ^ using 

« (42) 

Re{A,Bl) « -J!lk'[,„\-'cos{2k[r,i\), (43) 

we obtain the power spectrum at the late time 

If we assume rji cx /c~", the resulting power spectrum behaves P„ cx k'^~'^^. In order to fit to 
the observations, n ^ 1. 



5 

M 



For A;|?7i| < ^ 



1 

^^(^^^^•) ^ i^F^^°^(2A:|r/i|), (46) 



then the power spectrum is 



'P^k) ^ { ^X^k"'\^,r{2-cos{2k\r^^\)). (47) 



271 J 4if2 

Similarly for kn\r]i\ ^ M/H, if we assume ?7i oc Vh oc /c"^"*"^" and it is required to 
be n «i 1 in order for the spectrum to be scale invariant. Otherwise, it may not satisfy the 
observations. M/ H can also be constrained from the observations. 



IV. SPECTRUM OF THE RELIC GRAVITATIONAL WAVE 



In this section, we compute the number of gravitons produced at the de Sitter- radiation 
dominated phase transition in order to estimate the intensity of the gravitational wave due 
to the effect from the trans-Planckian physics for the detection at the gravitational wave 
detectors. 

With the scale factor during the radiation dominated phase {r]e < i] < rjeq) 

a{v) = ^iv-'^Ve), (48) 
where rje is the time when inflation ends, the mode solutions of are given by 

MV) = -^e-^'^ + 4^6^^=^ (49) 



and the coefficients are satisfied the normalization condition 

\ak\^-\(3k\^ = l. (50) 



Requiring the mode solutions ( 1491) during the radiation dominated phase and (1151) during 
the inflationary period to be continuous at 77 = r/e gives 

The spectrum of the relic gravitational waves is expressed in terms of the number of 
gravitons produced at the transition 10 1 



= '-P^ = ^iV. (53) 

Pcdlogf Pc 

where is the energy density of the stochastic background gravitational waves, / is the 
physical frequency and pc = 3H^Mp/{87i) is the critical energy density today. 
From (l52l) . N^^ is given by 



where we have used l^^p — = 1. Since kphys = ^ = 2'/r/, we can express k\rif>\ in terms 
of the physical frequency observed today / 



/ 

k\Ve\ = y, (55) 

Je 



where 



Here we have chosen the reference values for He as 10~^Mp and taken in order to reach to 
the numerical values in ( l56l) 

te = 1+ZegC^ 2.4 X lO^nohl, 

IHe 

teq ~ 4.1 X 10^°fio^/io^S 

Mp = 1.22 X lO^^GeV = 1.85 x lO^^Hz. 
Thus, the resulting spectrum of the relic gravitational wave can be written as 

He Vff 



hln,^ = 1.2n' X 10-1^ ( ) [jl (57) 



Similarly, A;|?7i| can be expressed by 



k\Vi\ = T' (58) 
Ji 



/i = .;r7^--^'f^V" (59) 



where 



27r{l + Zeq) \teq/ 

= e-^^^fe (60) 
and Afi = Hdt. If we assume Hi ~ He during inflation , since Afi ^ 1, we have always 

fl < fe. 

The modes having frequencies with f < fi exit the horizon during the trans-Planckian 
regime and the modes with fi < f < fe cross outside the horizon during the inflationary 
period in which the linear dispersion relations hold (I r/e I < \v\ < I'?!!)- Since the characteristic 
frequency corresponding to the radiation-matter equal time {rjeq) is feq — lO^^^^o^oHz, 
the modes that have a frequency / > feq today re-enter the horizon during the radiation 



dominated phase and the modes with / < f^q during the matter dominated phase. In this 
paper, we focus on the frequencies with / > feq. 

Since it is only meaningful to consider / <^ /e, we can approximated the produced number 
of particles as 



iV, 



1 
4 



\Bk\' - Re{AkBl 



(61) 



- 4 ■ 4' 

For the usual standard inflation scenario with an adiabatic initial vacuum (Ak = 1, Bk = 0), 
since A^^ oc f~^, the relic gravitational wave shows the fiat spectrum. Even if we consider 
the effect of the trans-Planckian regime, if rji is given by f|T7|) . as in fl22|) the relic spectrum 
of the gravitational wave is independent of /. 

Taking into account of the ignorance of the transition process around at rj = rji as in HL 
gravity (Sect. IIIip . in this section we will leave /i as a free parameter. 

If we define 



fc 



1/m 



M 



f: 



1) 



for / ^ /i, using the asymptotic form of Hi 



(2) 



f 5 



H' 

231). we obtain 



(62) 



\Bk\ 



Re{AuBl) 
Hence, the spectrum for f ^ fc is 



8\f 
1 



2m-\ 



- \ ^ \ COS 



/ 



2/ 

fi 



2m-\ 



(63) 
(64) 



h^o^g^ ~ 1.85 X 10 



-13 



10-Wp 



- + cos 



and for f <^ fc 



4.6 X 10" 



-14 



cos 



2/ 
/i 



2/ 

/i 



(65) 



(66) 



10-4Mp^ 

As a consequence, for / ^ /i, the shape of the spectrum changes at / = fc. If we assume 
/i oc /" as in Sect. Illlt the spectrum grows as H f > f^ and decreases as f-'^i^-'^) 

if f < fc- m = OT n = 1 gives a fiat spectrum. Note that the expression fl66|l is, in fact, 
not exact form, because, for /i ^ / ^ fc, z\ oc {f I fcT f I f\ and then the asymptotic form 
of the Hankel function ( l23l) is not appropriate for this range. The exact result can be found 
numerically, which is shown in Fig. [H 




FIG. 1: h^flgw versus / for different values of A/i, H/M and m 



For / <^ /i, the Hankel function is approximated as 



1 



Z Sm TTZ/ 



_T{l + iy)\2j r(l 



1 



if 2; < 1. 



Since -C 1, 2 term is dominant over z'^ term and then we have 



]\^\ 2(l+m) / 



f 



Re{AaBl) oc 



/ 



Using these coefficients, the spectrum becomes 



M 



2{l + m) / 

7 



(67) 

(68) 
(69) 

The spectrum for f <^ fi decreases y-9(i-")/2 which is independent of m, where we have 
assumed /i oc 

In Fig. [T], we plot the spectrum of the relic gravitational wave, h^Qgy^, with the physical 
frequency observed today for different values of A/i, H/M and m. m = 2 corresponds to 
HL gravity. The graph inflects at / = /c which is given in (162|) and also bends at / = /i. 
If we take M = 50 and H/M = 10-^ ~ lO^^Hz. As expected from ii f > fc, the 
spectrum increases with /"^ as / increases. And if /i < / < /c, it decreases for increasing /, 
but the slope is insensitive to the m as explained in below (l66l) . For / < /i, the spectrum 
decreases with /~^/^, independently of m, as / increases. 




FIG. 2: hlVtg^ versus /i with H/M = 10"^ for / = lOOHz. 



This figure seems to indicate the violation of the cosmological bounds from nucleosynthe- 
sis or CMB in the low frequency range 



101 ] or the present ground-base gravitational wave 
detector (LIGO, VIRGO etc.) in the range lOHz < / < lO^Hz. Even if m = 0, since 
h^flgw oc /~^^^ in the low frequency range, it might violate the cosmological bounds unless 
r]i is constrained by ( |T7I) . This implies that if /i oc /" where n should be determined from 
the transition mechanism, we can expect ^ 1 in order to satisfy the experimental and 
cosmological bounds. 

In Fig. [21 we plot h^Qg^^ as a function of /i at fixed / = lOOHz which is the frequency 
range of the ground based gravitational wave detectors. The spectrum decreases for increas- 
ing fi independently of m. For a given /i, as m increases, the amplitude of the spectrum 
also increases. As a result, /i could be constrained from the experiment, if the gravitational 
waves are detected. 



V. SUMMARY AND DISCUSSION 



We have estimated the spectrum of the relic gravitational wave if the modified dispersion 
relations in the trans-Planckian regime is given by the Corley-Jacobson dispersion relations. 
The Corley-Jacobson type dispersion relations could be obtained from the recently suggested 
quantum gravity model called Hofava-Lifshitz gravity. Hofava-Lifshitz gravity was proposed 



as a model to be renormalizable quantum gravity theory in the UV region and Einstein 
gravity is expected to be recovered at the IR hmit. m = 2 in the Corley-Jacobson dispersion 
relations can be considered as HL gravity. 

In obtaining the spectrum of the relic gravitational wave due to the effect of the trans- 
Planckian physics, HL gravity phase is assumed to be followed by the inflationary phase 
in Einstein gravity and how to occur the transition seems to play an important role. As 
in many literatures dealing with the trans-Planckian effect, if the instantaneous transition 
time is chosen by the cut-off momentum, the resulting spectrum in the late time is scale 
invariant and leads to the correction term by H/M. Since it is, however, not clear how 
the transition from HL gravity in the UV regime to the IR limit i.e. Einstein gravity takes 
place, we leave /i, the characteristic transition frequency, or r/i, the characteristic transition 
time, as a free parameter which is assumed to keep the information about the transition. 
We have found that in order to satisfy the experimental bounds and cosmological bounds 
from the nucleosynthesis or CMB, if we assume /i oc /" or rji oc n ~ 1. 

Although it is still controversial if HL gravity is a complete theory, it would be interesting 
to investigate whether HL gravity as like any other trans-Planckian physics can provide any 
detectable effects on the observations. And if HL gravity can provide the seeds for the 
large scale structure without inflation as well as solve the problems of the standard Big 
Bang cosmology, it would be also interesting to find some observational signature due to the 
tensor perturbations as well as the scalar perturbations in HL gravity without inflation [l^] . 
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